Decoherence suppression of open quantum systems through a strong coupling to 

non-Markovian reservoirs 
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In this paper, we provide a mechanism of decoherence suppression for open quantum systems in 
general, and that for " Schrodinger cat-like" state in particular, through the strong couplings to non- 
Markovian reservoirs. Different from the usual strategies of suppressing decoherence by decoupling 
the system from the environment in the literatures, here the decoherence suppression employs the 
strong back-reaction from non-Markovian reservoirs. The mechanism relies on the existence of the 
singularities (bound states) of the nonequilibrium retarded Green function which completely deter- 
mines the dissipation and decoherence dynamics of open systems. As an application, we examine 
the decoherence dynamics of a photonic crystal nanocavity that is coupled to a waveguide. The 
strong non-Markovian suppression of decoherence for the optical cat state is attained. 

PACS numbers: 03.67.Pp; 03.65.Yz; 42.50.Dv; 05.70.Ln 



I. INTRODUCTION 

How to protect quantum states away from decoherence 
is one of the most challenge topics in quantum informa- 
tion processing and modern quantum technology. During 
the past two decades, many schemes have been theoret- 
ically proposed and experimentally realized to suppress 
the decoherence in quantum information processing [J- 
6] . On the other hand, due to the significant development 
of the nanotechnology during the past decade, various 
quantum devices with high tunabilities, such as nanome- 
chanical oscillator or superconducting qubit strongly cou- 
pled to cavity 0, M , trapped atom coupled to an en- 
gineered reservoir [S | , arrays of coupled nanocavities in 
photonic crystal [10 1 etc., can be engineered. In these 
quantum devices, the strong coupling between the sys- 
tem and the structured reservoir and the resulting non- 
Markovian back-action play an important role in the ma- 
nipulations of quantum coherence. 

In this work, we shall provide a general mechanism 
of decoherence suppression for quantum systems coupled 
strongly to non-Markovian reservoirs. Contrary to the 
ordinary means of suppressing decoherence via dynami- 
cally decoupling of the system from the environment 0- 
[g], we employ the strong non-Markovian back- reaction 
from the environment to suppress the decoherence of 
quantum states. We show in general that when the non- 
Markovian back-reaction is strong enough, the decoher- 
ence of quantum states can be largely suppressed. In 
particular, we examine the time evolution of the Wigner 
function for a mesoscopic superposition of two coherent 
states, and demonstrate that the decoherence of such a 
mesoscopic superposition state can be suppressed due to 
the strong non-Markovian back-reaction from the envi- 



ronment. 



II. EXACT MASTER EQUATION 

The dynamics of open quantum systems are described 
by the reduced density matrix which can be obtained by 
tracing over all of the reservoir degrees of freedom form 
the total system p(t) = tr[p to t(t)], where ptot(^) is the 
total density matrix of the system plus its reservoir. The 
exact master equation of the reduced density matrix p(t) 
for an open system, such as a cavity in quantum optics, 
a defect (nanocavity) in photonic crystals or a quantum 
dot in nanostructures, etc., coupled to a general non- 
Markovian reservoir has been derived recently [TTl - fT3L I20I ] 



dp(t) 
dt 



-[H'(t), p(t)] + 7(t)[2ap(£)a t - p(t)a) a - a) ap(t)] 



- 7(t)[ap(t)a^ + a) p(t)a — ajap(t) — p(t)aa f ] 



(1) 



where H'(t) = uj'(t)a^a is the renormalized Hamiltonian 
of the system with the renormalized frequency uj'{t) — 
— Im[u(i)u _1 (t)]. The time-dependent coefficients jit) = 
-Refu^)^ 1 ^)] and j(t) = v(t) ~ 2v(t)Re[u(t) U - 1 (t)} 
incorporates all of the dissipations and fluctuations in- 
duced from the coupling to the reservoir. The function 
u(t) is the nonequilibrium retarded Green function of the 
system satisfying the following equation: 



ii(t) + iuj c u(t) + g(t- t)u(t) = 



(2) 



subjected to the initial condition u(£o) = 1> and the 
nonequilibrium thermal fluctuation is characterized by 
the function v(t) which is given by 
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V(t) 



dr dT'u*(Ti)g(n - t 2 )u(t 2 ). 

Jtn Jtn 



(3) 
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By introducing the spectral density of the reservoir 
J(uj) = 2jT^2 k \Vk\ 2 S(uj — ook) where is the coupling 
between the system and the reservoir, the time correla- 
tion functions g(r — r') and g(r — r') in Eqs. (J5J3]) are 
given by 



g(r-r')= r^j (w )e-M~0 



(4) 



duj 



9(t-t')= — JHn( W ,T)e— ') , (5) 



2tt 



which characterize all the non-Markovian back-reactions 
of the reservoir, and n(uj,T) = ehbJ /k B T — j" is the average 
particle number distribution in the reservoir at the initial 
time io- 

The decoherence dynamics of quantum states can be 
studied by examining the evolution of the correspond- 
ing Wigner function. With the help of the exact master 
equation ((T|), the exact Wigner function of an arbitrary 
quantum state at arbitrary time t in the complex space 
{z} is found: 



W(z,t) = J dn(a Q )d^{a / Q )(a \p{t )\a' }^(z,t\a 0l a' Q *,t ), 

(6) 

where |a) = e aa |0) is the coherent state, d/i(a) — 
d " da e ~ ' a ' is the integral measure of the Bergmann com- 
plex space, p(io) is the reduced density matrix of the ini- 
tial state, and the propagating function 3(z, t\ao, a'J , to) 
is given by 

Z(z,t \a ,a *,t Q ) = W Q {z,t) exp {z*tt(t)u{t)a 

+ a'*u*(t)n(t)z + a , *(l-\u(t)\ 2 n(t))a } , (7) 

where = and W (z,t) = exp [n{t)\z\ 2 ]. 

To concentrate on quantum decoherence, we examine 
the time evolution of a mesoscopic superposition of two 
coherent states moving in opposite directions, called as 
the "Schrodinger cat-like" state or the optical cat state 
in the literature 14]: \4>) = N(\a) + \ - a)), where N = 
cosh \a\ 2 is the normalization factor. As a result 
of Eq. (|6]), the time-evolution of the Wigner function for 
this cat state is given by 



W(z, t) = W a {z, t) + W~ a (z, t) + Wi [z, t) 



with 



W±M,t) = iv2^ e M 2 e fi(*)l^(tH 2 ; (8b ) 
Wj(z,t) = 2N 2 ^er^\{e-^~ u ^' n ^ z+u ^ a n . 

TT 

(8c) 

In Eq. ([5]), the first two terms are the Wigner functions 
for the initial coherent states \a) and | — a) respectively, 
the third term is the interference between them. The 



quantum coherence of the cat state can then be charac- 
terized by the fringe visibility function 



F(a,t) 



1 



Wi{z,t) | peak 



2 VWa(-M)|peakW / -a(M)lpeak 

K*)l 2 



exp 



{-2M 2 (i 



1 + 2v{t) 



>} 



(9) 



which ranges from unity to exp (— 2|a| 2 ) for full coher- 
ence to complete decoherence. As shown by Eqs. ([3]), 
(|8]| and ([9|), the nonequilibrium retarded Green function 
u(t) completely determine the dynamics of the quantum 
decoherence of the system. Equation ([2]) alone can also 
give the exact solution of atomic systems involving only 
single excitation (single photon process) at zero temper- 
ature [IE [l6| . 



III. GENERAL MECHANISM OF 
DECOHERENCE SUPPRESSION 

The solution of the retarded Green function can be 
obtained by the inverse Laplace transformation [l7Hl9l | 
u 0) = 2^1 Ib dsu{s)e st , where u{s) = ls _ u l_^ s) and the 
Bromwich path B is a line Re(s) = const > in the half 
plane of the analyticity of the transformation. The self- 
energy S(s) = J T^~L is the Laplace transformation 
of the correlation function (j4j. Consider a spectral den- 
sity ranged from w e to infinity, e.g. cj e = for Ohmic, 
super-Ohmic and sub-Ohmic reservoirs, etc. The self- 
energy is then not defined on the segment of the imag- 
inary axis s = —ilj with u) > w e , while s = —iui e is a 
branch point. Near the imaginary axis, the self-energy 
function can be separated into real and imaginary parts 
by the relation lim^o u ± iTI = V ^^iTi5(uj) so that E(s = 

-iu ± 0+) = A(w) T with A(w) = V f°° ¥-^, 
and V denotes the principal value. The analytic prop- 
erties of the transformed retarded Green function u(s) 
determine completely the decoherence dynamics of the 
system. 

In the very weak coupling regime, the self-energy func- 
tion is dominated near the pole s = —ilo c . The functions 
A(w) and J(u>) can be approximated by A c = A(w c ) 
and J c = J(uj c ). The resulting retarded Green func- 



3a) tion becomes u(t) = e~ 



2 * with the shifted frequency 



lj' c = uj c + A c . The retarded Green function experiences 
an exponential decay with the decay constant J c /2, which 
reproduces the Born-Markov result [HI HO]. Thus, u(t) 
will eventually be damped to zero and the fringe visibil- 
ity will be decayed to exp (— 2|a| 2 ), namely the quantum 
coherence is totally lost. 

However, as the coupling increases, the variation of the 
self-energy away from the pole — iuj c becomes significant, 
and the decoherence dynamics of the system is then to- 
tally different from the Born-Makov limit. Especially, 
there exists an isolated pole s = on the imaginary 
axis outside the branch cut, i.e. — lj c = A(57) with 
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FIG. 1: Integration contour of the inverse laplace transform of 
u(s). The Red line on the imaginary axis is the branch cut. 
—iQ. is the pure imaginary pole. The integration along the 
solid (dashed) curve is made on the first (second) Riemannian 
sheet. Si and s'i are the poles of u(s) on the first and second 
Riemannian sheet in the Re(s) < half-plane. 



fl < uj e , which leads to a dissipationless dynamics of the 
system. The exact solution of the retarded Green func- 
tion can be obtained by the inverse Laplace transform 
along the Bromwich path B as shown in Fig. [T] Since 
the closure crosses the branch cut Im(s) < — w e on the 
imaginary axis, the contour is necessary to pass into the 
second Riemannian sheet in the section of the half plane 
with Im(s) < — w e , where it remains in the first Rieman- 
nian sheet in the sections Im(s) > — w e in the half plane 
Re(s) < 0. To properly close the contour, it is neces- 
sary to turn around the branch point — iu e , following the 
Hankel paths h to enter and leave the second Rieman- 
nian sheet, as shown in Fig. Q] The exact propagating 
function can be obtained by means of the residue method 



u(t) =Ze- mt + Z l e-^-' l ^ t + Z 'it 



-'y.t—ibjA 



Q-iu 



— I ds[u n (s)-u\s)] 



(10) 



where Z is the residues of the bound state with the 
imaginary pole s = —iCl, Zi(Z'j) is the residues of the 
ith unstable states with the pole s = Sj = —7; — iuji 
(Sj = — 7! — iuj'j) on the first (second) Riemannian sheet 
which is the solution of is — oj c — E I,II (s) = with 
£ n (s) = E ! (s) + iJ(is). The last term is the contribu- 
tion from the contour along the Hankel path h (Fig. [T]) 
which is responsible to the nonexponential decay dynam- 
ics [13] ■ As a result [shown by Eq. (fTt)|)]. the retarded 
Green function shows the dissipationless dynamics due 
to the existence of the bound state. This means that the 
decoherence of the system can be suppressed through the 
strong non-Markovian coupling to a reservoir. The co- 
herence preservation in the cat state is also obvious by 
substituting Eq. (fT0|) into Eq. ([9]). It is straightforward 
to extend the above analysis to structured reservoirs with 
finite spectrum, as shown explicitly in the following dis- 
cussion. 
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FIG. 2: Graphical solutions of the imaginary poles of the 
retarded Green function, (a) Below the critical coupling, no 
solution outside the branch cut. (b) Over the critical coupling, 
two solutions outside the branch cut. 



IV. AN EXAMPLE FOR APPLICATIONS 

As an application, we apply the above general mecha- 
nism to the decoherence dynamics of a nanocavity (with 
frequency oj c ) coupled to a structured waveguide [its 
characteristic dispersion ujk — l<Jq — 2£o cos (&)]. The cou- 
pling strength between the nanocavity and the waveguide 

in photonic crystals is Vk = £ sin (nk) [HI, [lj| . The 

spectral density J(w) is then given by 



wo J 



\uj - w 1 < 2£ 
|w - w 1 > 2^0 



(11) 



where r\ = £/£o characterizes the strength of the cou- 
pling between the nanocavity and the structured reser- 
voir. From the above spectral density, the self-energy in 
ii(uj) can be exactly calculated 



S(s) = -v 2 [(a + iw ) - V(2£o) 2 + (s + iw ) 5 



(12) 



As the coupling strength exceeds the critical value 

i] c = 2 — ii££__f£°i j bound modes (the poles determined 

graphically in Fig. [2]) occur. As a result, when the cou- 
pling strength is below the critical coupling, no imaginary 
pole exists outside the branch cut, see Fig. [UJa). The so- 
lution of the retarded Green function shows a dissipative 
dynamics. However, when the coupling strength is larger 
than the critical coupling, one or two imaginary poles 
appear, see Fig. [2jb), and the solution of u(t) behaves 
dissipationless after a short time. 

To see explicitly the mechanism of decoherence sup- 
pression through the strong non-Markovian effect, we 
may look at the steady state solution of the nanocav- 
ity in the strong coupling regime. Consider the case the 
frequency of the nanocavity equals to the band center of 
the reservoir, i.e. w c = wq, the steady state solution of 
u(t) becomes 



(13) 



This shows that the retarded Green function is enveloped 
by the cosine function with the amplitude A(rf) 



and the frequency ui(j]) = 



£0 which corresponds to 




FIG. 3: (a) The exact numerical result of the retarded Green 
function u(t) and (b) the thermal fluctuation v(t) /n(ujo, T) in 
different coupling strength. The frequency of the nanocavity 
uj c is set to be the same as the band center of the waveguide 
coo- 



the energy exchange between the cavity and the reservoir. 
Fig. [3] shows the exact numerical result of the retarded 
Green function u(t) [see Fig. EJa)] and the normalized 
thermal fluctuation v(t)/n(ujQ,T) [i.e. Fig. Efb)] in dif- 
ferent coupling strength. Note that when the coupling 
rj > rj c = y/2, both the retarded Green function u(t) 
and the thermal fluctuation v(t) keep oscillating rather 
than damping. The oscillation indicates that the cavity 
keeps exchanging photons with the waveguide due to the 
strong non-Markovian back-reaction from the reservoir. 
The steady solution of the fringe visibility function of 
Eq. (|9|) at zero temperature simply becomes 



F(a, t) = exp {-2|a| 2 (l - A(rj) 



2 cos 2 



Hv)t}} ■ (14) 



Instead of full decoherence, the cat state keeps oscillat- 
ing in the strong coupling regime. The stronger the cou- 
pling strength is, the larger the degree of coherence can 
be maintained. Fig. 2] shows the periodic motion of the 
Wigner function for the cat state with the coupling r\ = 4 
and the temperature T = 0.5mK. As shown in Fig. 01 the 
interference of the cat state keeps oscillation in time. 

In fact, in the weak coupling regime, the fringe visibil- 
ity will eventually decay to e 2 1 I because of the deco- 
herence induced by the reservoir. Then all the coherence 
information of the cat state will be lost. At the same 
time, the larger the initial temperature of the reservoir is, 
the faster the decoherence processes. According to (|8bl) . 
as the the retarded Green function decays to zero, the 
two peaks of the Wigner function gradually spiral to the 
origin (see the movie for this Markovian time evolution 
given in [2l[) and the thermal fluctuation v(t) saturates 
to the equilibrium value n(uio,T) due to energy relax- 
ation. The cat state finally decays to a thermal state 
with the Wigner function 



W(z,t -> oo) 



7r[l + n(u> , T) 



exp[ 



2UI 



FIG. 4: Time evolution of the Wigner function of the meso- 
scopic superposition state. The coupling strength 77 = 4 and 
the cavity frequency equals to the band center of the waveg- 
uide uj c = (jJq, To = 2i{/(jJq. A movie for the above non- 
Markovian time evolution is given in [21] ■ 



In contrast, as shown in the previous analysis, when 
the coupling strength exceeds the critical coupling ?y c , the 
decoherence dynamics of the cavity field is totally sup- 
pressed. The fringe visibility, after a short time decay, 
oscillates above the value of e _2 ' Q ' for all the time. In 
other words, the coherence of the cat state goes to dead 
and birth repeatedly. In addition, according to Eq. (|8bj) 
and the stationary solution of Eq. (|13|) in the strong 
coupling regime, the two peaks of the Wigner function 
would keeps spiralling in and out of the origin with the 
=£0 due to the energy exchange between 



1 + n(u ,T) 



(15) 



frequency —fl= 

V'' 2_1 

the system and the reservoir, see the movie for this non- 
Markovian time evolution in [21]. Thus, the cavity field 
would never be thermalized by the reservoir and the de- 
coherence of the system is significantly suppressed. 



V. CONCLUSION 

In conclusion, we have shown through the exact master 
equation that the nonequilibrium retarded Green func- 
tion can completely determine decoherence dynamics. 
From analytic properties of the retarded Green function, 
we provided a general mechanism of decoherence suppres- 
sion through the strong non-Markovian back-reaction 
from environments. In particular, when the coupling be- 
tween the system and the reservoir exceeds a critical cou- 
pling, the bounded modes (the imaginary poles of the 
retarded Green function) leads to a dissipationless dy- 
namics such that decoherence can be largely suppressed, 
as a strong non-Markovian memory effect. This generic 
behavior is explicitly demonstrated through the decoher- 
ence dynamics of the cat state. Since the nonequilibrium 
retarded Green function is well-defined for arbitrary open 
quantum system, the mechanism presented in this work 
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